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ABSTRACT 


Singularly  perturbed  second  order  elliptic  partial  differential  equations 
with  Neumann  boundary  conditions  arise  In  many  areas  of  application.  These 

problems  rarely  have  smooth  limit  solutions.  In  this  paper,  w»  characterize 

r 
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the  limit  solution  for  a  wide  class  of  such  problems.  We-  also  give 
abstract  rate  of  convergence  theorem  and  apply  the  abstract  theorem 
finite  difference  approximations. 
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§1.  INTRODUCTION. 

In  this  pap«r,  we  study  the  singular  perturbation  problem  for  partial 
differential  equations  «dilch  have  the  form: 

-eAUg  +  H(x,Ug,7Ug)  “0,  x  e  R  , 

3u. 


(NP^) 


Y{x),  X  e  3R  , 

where  R  Is  a  bounded  domain  in  n  is  R's  outward  unit  normal,  u^ 

Is  a  scalar  unknown  and  H  Is  a  continuous  function  on  R  x  r  x  r*^.  one 
ai^llcatlon  that  motivates  the  study  of  singular  perturbation  problems  of  the 
form  (NP^)  Is  found  In  the  theory  of  optimal  stochastic  control.  There,  H 
depends  on  the  deterministic  part  of  a  stochastic  ODE,  a  control  space  and  a 
specified  cost  function,  u^  can  be  Identified  as  the  optimal  cost 
function.  The  positive  parameter  e  In  (MP^)  can  be  regarded  as  the 
Intensity  of  noise  in  the  dynamics  equation.  Control  problems  whose 
trajectories  reflect  at  a  boundary,  give  rise  to  Neumann  problems  of  the  type 
studied  here;  see  [1]  or  [17]  for  a  detailed  treatment  of  this  topic.  One 
could  ask  for  Instance,  Is  the  optimal  cost  function  of  a  stochastic  control 
problem  related  to  the  optimal  cost  of  its  associated  deterministic  problem? 
Are  the  two  close  In  any  way  when  the  noise  Is  small? 

As  e  'f  0,  it  Is  well  known  that  solutions  of  do  not  generally 

converge  to  a  classical  solution  of: 
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This  exan^le  Is  clearly  overdetermlned,  and  here  the  data  at  0  is  not 
compatible  with  the  data  at  1.  For  this  reason,  a  more  general  class  of 
solutions  to  (NPq)  must  be  sought* 

A  new  notion  of  continuous  weak  solutions  to  equations  of  Hamilton- Jacobi 
type  has  recently  been  Introduced.  In  [4]  and  [5],  M.  G.  Crandall  and  P.  L. 
Lions  have  developed  techniques  that  have  been  extremely  successful  In 
establishing  a  number  of  new  results  concerning  continuous,  but  not 
necessarily  differentiable,  weak  solutions  to  first  order,  fully  nonlinear, 
partial  differential  equations.  In  their  work,  Crandall  and  Lions  have 
utilized  the  "vanishing  viscosity  method",  so  neuaed  because  of  the  link  to  the 
classical  technique  of  vanishing  viscosity  from  fluid  mechanics,  and  they  show 
that  the  method  of  vanishing  viscosity  gives  rise  to  a  specific  notion  of  a 
"viscosity"  weak  solution. 

In  [15],  and  here  as  well,  the  notion  of  a  viscosity  solution  for  the 
generally  over  determined  Neumann  problem  (NPq)  is  given  and  Is  shown  to 
Include  all  L  c-llmlts  of  solutions  to  (NP^).  All  L  c**llmlts  are  shown 
to  satisfy  the  so-called  viscosity  Inequalities  of  Section  2;  additional 
details  In  this  direction  can  be  found  In  [15].  Remarkably,  (with  additional 
hypotheses  of  course),  these  viscosity  Inequalities  uniquely  determine  all 
such  limits.  In  Section  3,  we  Introduce  what  we  call  "approximate  viscosity 

Inequalities"  and  we  show  there  that  any  reasonable  approximate  viscosity 
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solution  is  approximately  equal  to  the  viscosity  solution  of  (NPg).  More 
precisely,  we  give  an  abstract  minimal  rate-of -convergence  theorem.  Theorem  2, 
for  approximate  viscosity  solutions  to  (MPq).  We  also  show  that  this  rate  is 
essentially  sharp.  A  particular  application  of  Theorem  2  gives  an  easy  to 
determine  measure  of  how  far  the  solution  of  (NP^)  can  be  away  from  the 
viscosity  solution  of  (NPq).  In  Section  4,  the  abstract  rate-of-convergence 
theorem  of  Section  3  is  applied  to  numerical  approximations  which  are  obtained 
from  a  class  of  finite  difference  schemes.  Moreover,  we  show  in  Section  4 
that  these  schemes  have  "computable”  solutions  and  we  motivate  how  they  can  be 
obtained. 

The  reader  is  encouraged  to  see  [2],  [18]  and  [19]  where  similar  results 
as  those  above  are  obtained  for  divergence  form  singular  perturbation  problems 
with  mixed  or  Dirichlet  boundary  conditions.  Also  see  [6],  [7]  and  [21]  for  a 
further  treatment  of  approximations  for  time  dependent  Hamilton-Jacobi 
equations  without  spatial  boundaries. 

52  VISCOSITY  LIMIT  SOLUTIONS. 

As  mentioned  in  the  previous  section,  as  e  4’  0,  the  corresponding 

solutions  to  (I4P^)  do  not  in  general  converge  to  a  classical  solution  of 

(NPq).  In  this  section  we  offer  a  characterization  of  viscous  limits  to  (NPq) 

and  we  show  that  this  characterization  often  allows  for  only  one  solution  in 

the  class  of  continuous  functions.  Throughout,  we  shall  assume  that  0  is  a 

bounded  domain  in  which  has  a  boundary  30.  Ihe  outward  normal  of 

0  at  a  point  x  e  30  will  be  denoted  by  n(x)  and  we  write  the  outward 

normal  derivative  of  <p  at  x  8  30  as  (x). 

dn 

We  should  like  to  mention  that  previous  to  the  writing  of  this  paper 
P.  L.  Lions  had  introduced  the  same  viscosity  characterization  of  solutions  to 
(NPO)  we  give  below;  see  [IS].  For  this  reason,  we  borrow  much  of  the 


notation  and  hypotheses  of  [15]  and  In  this  section  we  osiit  all  proofs  but 
those  which  motivate  the  results  of  the  next  sections. 


We  now  state  the  vicosity  characterization  (see  Proposition  1)  of 
continuous  weak  solutions  to  (NPq). 

Definition  1 .  Suppose  that  H(XfU«p)  €  C(fl  x  R  x  s^)  and  u(x)  6  C(()).  ffe 
say  that: 

(a)  u(x)  la  a  viscosity  subsolutlon  of  (MPq)  If  for  all  test  functions 

If  e  C^(R^)  with  >  y(x),  we  have 

H(Xg,u{xQ),7q>{XQ))  <  0  , 

where  Xq  e  n  satisfies 

u(Xq)  -  <p(Xq)  ■  max  (u(x)  -  sCx))  . 
xen 

(b)  u(x)  Is  a  viscosity  supersolution  of  (HPq)  If  for  all  test 

functions  <f  6  with  <  Y(x)»  %»o  have 

H(Xg,U(Xg),7s(XQ))  >  0  , 

where  Xq  e  ^  satisfies 

u(x-)  -  ifiXf.)  -  min  (u(x)  -  s(x))  . 

xen 

(c)  u(x)  Is  a  viscosity  solution  of  (NPq)  If  It  satisfies  both  (a)  and 
(b)  above. 

The  fact  that  our  test  functions  are  required  to  satisfy  ^  Y(x), 

(resp.  ^  y[x)),  In  our  definition  of  viscosity  subsolution,  (resp. 

supersolution),  may  at  first  seem  superfluous.  Ihls  Is  however,  precisely  the 
mechanism  that  "sees"  the  Neumann  boundary  conditions  when  vanishing  viscosity 
Is  taken  into  account;  (see  Proposition  1  below). 

Remark  2. 1.  Any  solution  of  (NPq)  Is  also  a  viscosity  solution.  This 

fact  Is  nontrivial  only  for  the  case  when  max(u  -  s)  or  mln(u  ~  s)  Is 
attained  for  some  xO  6  3fi.  To  see  that  u  must  indeed  be  a  viscosity 


subsolution,  take  an  arbitrary  e  6  C^(R^)  with  ^  Y(x).  First  choose 

a  sequence  (e.)/  such  that  for  every  m,  o_(Xn)  •  e(Xn),  o_(x)  >  e(x)  for 
B  m  u  u  m 

d« 

X  V  xnf  T-  (x)  >  -r^Cx)  and  with  a  a  In  C'  as  m  •»■  ».  We  then  have 
w  3n  3n  ra 

that  (u  -  e)(xQ)  >  niax(u  *•  and  Xq  is  the  point  where  the  strict 

maximum  of  u  -  e  Is  attained.  Next,  for  any  fixed  m,  choose  a  sequence 

n 

de’^  de 

such  that  ^  with  e”  in  as  n  +  o». 

n  dn  dn  n  n 

Denoting  by  x^^  the  points  where  max(u  **  Is  attained,  we  must  have 

that  x^  Xq  as  n  «.  This  Is  true  because  (u  -  ^  strict 

d  n 

maximum  of  u  -  e  .  For  x  6  3fl,  we  also  have  that  -^(u  -  <  0, 

m  dll  a 

which  Implies  x^  6  (2^.  Therefore,  since  now  x^  Is  an  Interior  maximum  of 

u  -  e",  Vu{x  )  •  Ve'^tx  ),  and  so  by  taking  limits  we  have 
n  n  n  n 

H(x  ,u(x-) ,Ve(x  ))  ■  lira  H(x  ,u(x  ) ,7e”(x  ) )  -  0  . 

0  0  0  n  n  m  n 

ra,n 

Remark  2.2.  Obviously,  the  converse  of  Remark  2.1  Is  false.  That  Is,  a 

smooth  viscosity  solution  need  not  satisfy  the  boundary  conditions  of  (NPq). 

2  - 

Proposition  1.  Let  u^  6  C  (fi)  be  a  solution  of  (NPg)  and  suppose  that 
2  - 

e  e  C  (fl) .  Then: 

(a)  For  ^  T(x)  and 

that  H(xQ,u^(XQ),7e(xQ))  <  zLfp{x^), 

(b)  For  |^(*)  <  Y(x)  and  u^(Xq) 
that  H(xQ,Up(xQ),7e(xQ))  >  eAe(Xg). 

If  In  addition,  we  have  that  u^  >  u  In 
then: 

(c)  u  «  11m  u  Is  a  viscosity  solution,  that  Is,  u  satisfies 

e 

Definition  1c. 

The  proof  of  Proposition  1  can  be  found  in  [15],  however  the  Interested 
reader  can  easily  reproduce  Its  proof  by  taking  limits  as  in  Remark  2.1. 


-  «(Xn>  -  max  (u„(x)  -  (p(x))  we  have 
xen 


-  i))(x-)  -  min  (u  (x)  -  (p(x))  we  have 

xejj  ^ 


0»  M 

L  (n)  for  some  sequence  e  'V  0, 


Before  stating  the  main  result  of  this  section/  we  give  a  simple  lemma* 
Lemma  1 .  Ijet  Q  be  a  bounded  domain  in  having  a  boundary  3(2. 

Then: 

(a)  niere  exists  a  constant  <  <"  such  that  for  all  x  €  3(2, 


Proof t  =  w(x)  =0  would  suffice  In  the  case  of  convex  (2.  For  nonconvex 

(2,  (a)  Is  shown  In  [12].  (b)  can  be  shown  by  constructing  a  particular 

example.  Under  the  hypotheses  of  the  lemma.  It  Is  known  that  the  distance 
function  d(x)  =  d(x;3(2)  Is  In  a  neighborhood  of  3(2;  [23],  [10].  Tha 

Is,  d(x)  e  C  ((2^),  where  (2^  *  {x  6  (2  ;  d(x)  <  x}  and  x  >  0  Is  chosen 
sufficiently  small.  Set  0  <  Xg  <  x  and  verify  that 


Now,  consider  the  following  set  of  assumptions. 

Assumption  A:  H(x,u,p)  Is  strictly  increasing  In  u  for  all  x  e  (2  and 
uniformly  for  p  6  Rd.  That  is,  for  all  R  >  0  and  -R  <  v  <  u  <  R,  there 


exists  a  >  0  such  that 


H(x,u,p)  -  H{x,v,p)  >  p  (u  -  v)  . 


Assumption  B;  Let  a,0  6  R*^  satisfy  |a|,|3|  <  max(C^/C)>  where  is  as 

defined  In  the  previous  lemma*  Then,  for  all  such  a, 6  and  all  x  e 
X  ^  e  all  |u|  <  R  and  any  X  >  1,  assume 

1h(x  +  C,u,XC  +  j  Ul^a  +  0(5))  -  H(x,u,X5  +  |  |5|^e)|  <  <*)j^(X|5l^  +  |5|)  , 

where  u)  (s)  Is  some  function  such  that  11m  u)-(s)  =<  0. 

S'i'O 

Remark  2.3.  Assumptions  A  and  B  are  standard;  see  [5],  [8]  and  [15].  In 
the  following  theorem.  Assumption  B  may  always  be  relaxed  so  that  a  =  3  =  0 
and  0(5)  =  0  except  for  x  In  a  neighborhood  of  30.  Assuming  additional 
regularity  on  the  class  of  solutions  allows  Assumption  B  to  be  neglected 
entirely. 

Theorem  1 .  Suppose  that  H(x,u,p)  6  0(5  x  R  x  R^)  and  that  it  satisfies 
Assumption  A  above.  Let  u  6  0(5)  be  a  viscosity  subsolution  of  (NPq)  and 
let  V  e  0(5)  be  a  viscosity  supersolution  of  (MPq).  Finally,  assume  one 
from  the  following  three  sets  of  hypotheses: 

(1)  n  Is  convex  and  Assumption  B  holds  with  a  =  3  =  0. 

(11)  Assumption  B  Is  satisfied. 

(Ill)  Either  u  or  v  is  Lipschltz  continuous. 

We  then  have  that 

max  (u(x)  -  v(x) )  <  0  . 

xen 

Obviously,  establishing  this  result  would  imply  that  a  viscosity  solution 
to  (NPq)  Is  unique  In  the  specified  class  of  functions. 


Proof : 

Given 

a  6  >  0, 

define  the 

function 

i|>^(x,y)  by 

(2.1) 

<t'*^(x,y)  * 

p(x)p(y) |x 

-  y|^/i5  t 

where 

P(x)  = 

exp(w(x) ) 

and  w(x) 

satisfies 

the  second  conclusion  of 

IJemma  1.  For  x  6  dn  and  emy  fixed  yg  e  u,  observe  that 
3  .6,  .  .6,  .r3w.  .  .  '  yo>*", 

^  4.  (x,yj,)  .  (x.yQ){^(x)  +-- - , 

X  I*  "  yfll 

and  Lenuna  1  implies  that  the  bracketed  term  above  is  nonnegative.  Therefore, 

h  ^  ° ' 

X 

and  similarly,  for  y  e  and  any  fixed  Xq  e  H 

Ik-  ^  °  • 

y 

Now,  choose  il*  ®  C^(S)  such  that  (x)  *  Y(x).  By  the  construction  eUaove, 
we  have  that  for  any  fixed  yQ  6  J5 

(|»^(x)  =  +  4(<x)  , 

is  an  admissible  test  function  according  to  Definition  la  and  similarly  for 


fixed  xq  e  n 


t,^(y)  o:  _,j)  (XQ,y)  +  It>(y)  , 


is  admissible  according  to  Definition  1b. 

The  next  step  is  to  note  the  obvious  inequality 

(2.2)  max  (u(x)  -  v(x))  <  max  (u(x)  -  v(y)  -  (♦^(x,y)  +  i|»(x)  -  ♦(y)))  • 

xeii  xen 

yen 

We  denote  by  x^,  y^  the  points  in  n  where  the  right  hand  side  of  (2.2)  is 
attained  and  we  rewrite  (2.2)  as 

(2.3)  max  (u(x)  -  v(x) )  <  u(Xg)  -  v(y,)  -  (♦  (Xg,yg)  +  <l'(Xg)  -  . 

xen 

Using  (2.3),  we  easily  arrive  at 

't>*^(xg,yg)  <  |u(xg)  -  u(yj)|  +  ('('(Xg)  -  '('(yg)!  . 


(2.4) 


’(Xx,y,)  <  |v(x^)  -  v(y,)|  +  |4'(x.)  -  4'(y;t)|  » 


and  recalling  the  definition  of  ^  (x,y),  (2>4)  gives  us  that 


(2.5)  ^  const. /6  . 

Furthermore,  since  u,  (or  v),  and  <|»  are  continuous,  (2.4)  combined  with 

(2.5)  shows  that 

(2.6)  11m  (t^(Xg,yg)  =  0  . 

6'(’0 

The  object  now  is  to  show  that  the  right  hand  side  of  (2.2)  can  be  made 
arbitrarily  small  by  choosing  6  sufficiently  small.  From  above,  we  see  that 
the  test  functions  defined  as 

'^^(x)  -  v(y^)  -  'Ky^)  +  i|i^(x,yj)  +  iKx)  , 

“  <t'*^(Xg,y)  +  iKy)  , 

are  admissible  according  to  Definition  la  and  Definition  1b  respectively. 
Inserting  these  into  Definition  1,  and  using  the  fact  that  u  is  a  viscosity 
subsolution  and  v  is  a  viscosity  supersolution ,  allows  us  to  conclude  that 


H(xg,u(x^),V^()>^(xg) )  <  0  , 

and 


H(y5,v(yg),Vy<|>2(y5))  >  0  , 


because  x^  satisfies 

u(xg)  -  4>^(x^)  *  max  (u(x)  -  i))^(x))  , 

xeH 

and  y^  satisfies 


v(yj)  “  't>2<y5^  “  "‘is  “  4'2^y))  • 

yen 

Combining  the  inequalities  above  and  rearranging,  we  obtain 


(2.7) 


H(Xg,u(xg),Vj^^^(Xg))  -  H(x5,v(yg),V^(t>^(Xg)) 

<  H(yg,v(yg),Vy4>2(y5))  -  H(  Xg  ,v(yg ) ,  ^  (Xg  ) )  . 


By  a  direct  calculation,  the  right  hand  side  of  (2.7)  can  be  written  as 


(2.8) 


where 


H(y^,v(y^),X(x^  ■  T  1*6  "  ^5!^^  ViKy^)) 

-  H(Xg,v(y^),A(x^  '  ^6^  1  1*6  "  ^5!^“  Vi|/(Xg))  , 


X  =  2p(xg )p(y^)/6 
a  =  Vw(Xg) 

B  =  Vw(yg)  . 


(Recall  from  Lemma 
p(x)  =  p(y)  =  1 ) . 

(2.9) 


1  that  if  (2  is  convex  we  may  assume  that  a  =  B 
Assumption  B  allows  (2.8)  to  be  bounded  above  by 


0  and 


where  R  =  max( lu|„, | v|„) . 

To  con^lete  the  proof  of  conclusions  (i)  and  (ii),  we  again  use 
inequality  (2.3)  to  write 


max  (u(x)  -  v(x))  <  u(x^)  -  v(y^)  +  U(Xg)  -  <|'(yx)|  , 

xeo 

which  is  bounded  above  by 

(2.10)  max((u(xg)  -  v(yg)),0)  +  |iKxg)  -  iKy^)!  • 

Assumption  A  applied  to  the  left  hand  side  of  (2.7)  combined  with  (2.8)  and 
(2.9),  allows  us  to  bound  (2.10)  by 


(2.11)  -J-  •‘'jjCMxg  -  yj^  +  Ixg  -  y^l)  +  It(xg)  -  t(y5)|  • 

R 

2  5 

Recalling  that  ~  ^6  ^  =2^  ^*6'y^6^'  ^2.6)  along  with  (2.5)  show  that 

(2.11)  tends  to  zero  as  6  tends  to  zero;  therby  proving 


max  (u( x)  -  v(x) )  <  0  . 


To  establish  (ili),  observe  that  If  U/  (or  v) ,  Is  Lipschltz 
continuous.  Inequality  (2.4)  leads  to  an  improvement  of  (2.5).  That  is,  we 
may  conclude  that 

(2.12)  jxg  "  y^ 1  <  const  6  . 

This  lnq>roved  estimate  implies  that  the  p  term  of  H(*,*,p)  in  (2.8) 
remains  bounded.  Therefore,  conclusion  (lii)  follows  by  noting  the  uniform 
continuity  of  H(x,u,p)  on  a  compact  subset  of  H  x  r  x  . 

Remark  2.4.  The  combined  results  of  Proposition  1  and  Theorem  1  imply  that 

m 

the  family  of  solutions  to  (NP^)  is  relatively  compact  in  L  i  then 

«D 

lim  u  exists  in  L  .  For  further  results  concerning  the  compactness  of 
e+0 

{Ug}g>o,  see  [13]  or  [14]. 

§3.  VISCOSITY  APPROXIMATIONS  AMD  A  RATE  OF  CONVERGENCE. 

In  this  section  we  consider  the  rate  at  which  certain  approximations 
converge  to  the  viscosity  solution  of  (NPq).  We  show  in  a  precise  sense 
below,  that  if  an  approximation  "almost"  satisfies  the  viscosity  inequalities 
of  Definlton  1  then  the  approximation  is  "almost”  equal  to  its  associated 
viscosity  limit  solution.  The  abstract  rate  of  convergence  theorem  given  in 
this  section  is  then  applied  in  Section  4  to  particular  approximations 
generated  by  a  class  of  numerical  schemes. 

Before  making  a  precise  statement  of  "almost  satisfies  the  viscosity 
inequalities",  recall  the  definitions  of  the  test  functions  used  in  the  proof 
of  Theorem  1 : 

(3.1)  <|>^(x,y)  =  p(x)p(y)lx  -  , 

where  6  >  0,  p(x)  exp(w(x))  and  w(x)  satisfies  the  second  conclusion  of 
Lemma  1.  Also  recall  the  function  i|/,  which  satisfies 


(3.2)  1(/(X)  e  C^lS) 

•|^(x)  -  Y(x)  on  30  , 
and  the  specific  test  functions 

(3.3a)  ♦^(x)  “  ♦^(x,yjj)  +  ^(x)  , 

(3.3b)  <p2<yi  “  -'<'^(XQ,y)  +  tl>(y)  , 

where  xq,  yQ  are  arbitrary  fixed  points  in  S. 

We  now  give: 

Definition  2.  Suppose  that  H(x,u,p)  6  C(S  x  R  x  R^)  and  ®  C(5).  We 

say  that : 

(a)  u^  is  an  approximate  viscosity  subsolution  to  (NPq)  if  for  all  test 
functions  vdiich  have  the  particular  form  (3.3a),  we  have  that 

H(Xo,Ug(Xo),Vx<jti(Xo))  <  eAx(>,(XQ)  +  Ce  , 

where  xg  €  H  satisfies 

Ug(XQ)  -  't'lCxg)  =  max  (Ug(x)  -  ♦^(x))  , 

xeo 


and  C_  is  some  fixed  consteuit. 

(b)  Ug  is  an  approximate  viscosity  supersolution  to  (MPq)  if  for  all 
test  functions  <^2'  which  have  the  particular  form  (3.3b),  we  liave  that 

H(yo/Ug(yo),  Vy(^2(yo))  >  > 

where  yg  ®  ^  satisfies 

Ug(yo)  -  <l'2<Yg)  “  "^S  • 

xeO 

(c)  Ug,  is  an  approximate  viscosity  solution  if  it  satisfies  (a)  and  (b) 


above . 


9^^ 

In  the  proof  of  Iheoren  1  we  showed  that  (x)  >  Y(x)  and 

3*2  * 

(y)  <  Y(y)«  therefore  the  statement  of  Proposition  1  implies  that  if 

y  2 

is  a  C*  solution  of  (NP  ),  then  it  is  also  an  approximate  viscosity  solution 
of  (NPn)  as  defined  above* 


With  Definition  2,  we  now  state t 

Theorem  2.  In  addition  to  Assumption  A  of  Theorem  1,  assume  for  ease  of 
presentation  that  *(x)  =  0.  Furthermore/  assume  that  (NPq)  admits  a 
Lipschitz  continuous  viscosity  solution  u#  with  say  Lipschltz  constant  L, 
and  assume  that  H(x,U/p)  is  locally  Llpschitz  continuous*  Then,  for  any 
approximate  viscosity  solution  to  (NPq),  say  u^,  we  have  that 

|Ug  -  u|.  <  -1-  [(SdLjjLd  +  2CjjL)e)^/^  +  0(e)] 


where 


Jq  -  sup  |u  I.  , 
1>e>0 


i*  -  sup  _ 
x^ 

l«l% 


|H(x^,u,P^)  -  H{X2/U,P2)| 

^  |x,  -  x^l  +  -  P2I  ^  ’ 


The  definition  of  is  given  in  Lemma  1. 


Remarlt  3* 1  *  We  may  replace  the  assumption  that  (NPq)  admits  a  Llpschitz 
continuous  viscosity  solution  by  the  assumption  that  u^  is  Llpschitz 
continuous,  uniformly  in  e  >  0* 

Remar)c  3*2.  The  Interested  reader  can  easily  modify  the  following  proof  to 
include  inhomogenious  boundary  data  to  obtain  the  same  /c  rate  of 


i.w  pyu  ".i  "ji  \w  rm  W  rMVfurJ.v.  WJi « 'VTM'VJvm.’vinKTVTrimrrj^.'rjw^jw^ 


convergence.  Relaxing  the  hypothesis  on  H(x,u,p)  and  the  regularity  of  u 
can  also  be  done  to  obtain  a  more  general*  (and  slower)*  rate  of 
convergence.  This  however*  will  not  be  done  here. 

Proof  of  Theorem  2;  Mlmlc)clng  the  proof  of  Theorem  1*  we  arrive  that  the 
analogue  of  Inequality  (2.7): 


(3.4)  H(Xj*Ug(Xj)*V^^^(Xg))  -  H(Xj*u(yj)*V^+^(Xj)) 


<  H(yg*u(yg)*Vy(j»2(y5) )  -  H(X5  *u(y5 )  (Xg ) )  +  eAj^(fr^(Xg)  +  Oe  * 

where  ®nd  (3.3)*  (with  \)»  H  0),  and  Xg  emd  y^ 

are  as  In  (2.3).  Recalling  (2.4)  and  using  the  fact  that  u*  (or 
Iilpschltz  continuous*  we  have  that 

(3.5)  'l>^(x^,y^)  <  |u(x^)  -  u(y^)|  <  Ljx^  -  y^  |  * 

or 

which*  with  the  definition  of  t^*  gives  us  that 


(3.6) 


1*6  " 


P(Xj)p(y^) 


Furthermore*  a  direct  calculation  shows  that 

^  2Cjj^^x,y)  , 

where  Is  as  In  Lemma  1.  This  Inequality*  along  with  (3.5)  shows  that 

<3-7)  “  ^y^2‘y6>l  ^  ’  ^6  ■  * 

Returning  now  to  Inequality  (3.4),  we  use  (3.7),  Assumption  A  and  the 
fact  that  H  Is  (locally)  Llpschltz  continuous  to  obtain 


Wj^  max(  (Ug(xg)  -  u(yg)),0) 


<  Ly(i  +  2CjjL)*|xg  -  yal  +  • 


Calculating  Inserting  (3 >6)  Into  the  right  hand  side  above/  we  find 


(3.8)  max((Ug(xg)  -  u(yj)),0) 

0 

^  ,P(X5)P(yfi)  -/6m 

<  [2de( - 3 - )  +  )]  +  const. (e  +  e6)  , 

6  **  6 


M-;- 


*'*w'*J> 


' 


where  L  ■  Ljj(  1  +  ^C^L) .  By  setting 


A _ r2d  ,>11/2 


y-  (^e)V2  , 

f  T  T 


p<*6)p<y6)  LL 

A 

which  can  be  done  for  LL  0  by  the  continuity  of  the  left  hand  side  with 
respect  to  6,  we  minimize  the  bracketed  term  in  (3.8).  This  yields 

Up  max((Ug(x3)  -  u(y3))/0)  <  (8dLLe)^^^  +  const. (e  +  , 


and  using  the  fact  that 

max  (Ug(x)  -  u(x))  <  -  uCy^)  » 

xe(i 

as  done  in  the  proof  of  Theorem  1/  we  have  established  the  desired  result  for 
max(Ug  -  u) . 

An  identical  estimate  can  be  obtained  for  max(u  ^  by  a  similar 

argument  and  so  the  proof  of  Theorem  2  is  complete. 


Remark  3.3.  When  the  domain  tl  is  convex  and  \(»  =  0,  the  term  0(e)  in  the 
estimate  of  Theorem  2  is  precisely  Ce.  in  addition,  if  the  approximate 
viscosity  solution  is  the  solution  of  (NP^),  the  constant  C  is  zero. 

Next,  we  show  that  the  order  of  the  rate  of  convergence  obtained  above 
can  not  in  general  be  improved.  To  see  this,  consider  the  example 


T*  "e  ”  ° 


(3.9) 


'•Vv 

V  ^  ■ 

••  .<*« 


-(0)  -  0, 


1)  •  1 


The  exact  solution  of  (3.9)  is  given  by 


Ug(x)  =  /e  co8h(x//e)/sinh( 1//e) 

and  its  an  easy  exercise  to  show  that  u^  0  uniformly  as  e  0.  In  fact. 


one  easily  finds  that 


lu  -  0|^  *  /T  {l  +  0(exp(-2//e) )}  , 


which  is  exactly  the  order  obtained  Theorem  2.  We  should  mention  however, 
that  the  rate  constant  of  Theorem  2  is  not  the  best  possible. 

We  conclude  this  section  by  £malyzing  the  specific  example: 

d^u_  ,  du_  2 

-e  — — -  +  (—  - — )  +  u  *  0  , 

,2  '■2  dx  e 

dx 


(3.10) 


du-  du_ 


Setting  e  =  0  and  solving  the  reduced  differential  equation,  we  find  that 

2 

u  =  lim  u^  should  be  built  from  functions  having  the  form  -(x  -  c)  and 

c 

0.  The  objective  now  is  to  piece  things  together  in  a  way  such  that  the 
constructed  function  satisfies  the  viscosity  inequalities  of  Definition  1. 
We  have  three  basic  cases,  (which  depend  on  Tn)*  Set 


•  o  •  ^ 


Ut  (X)  =  -(  x  -  •=-)  , 


f  0  ■»2 

Ur(x)  =  -[x  +  - -  , 


and  note  that  u^  satisfies  the  left  boundary  condition  of  (3.10)  and  Ujj 
satisfies  the  right  boundary  conditions. 


Case  1 .  For  1  >  Yq  ^  0,  consider  the  candidate  limit  solution: 

Uj^(x)  0  <  X  <  yi 

Ui(x)  -  j  0  _<  X  <  1  -  y- 

^0 

Uj^(x)  1  -  y  <  X  <  1  . 

The  analysis  of  this  case  is  trivial  since  u^Cx)  is  a  classical 
solution  to  the  reduced  problem.  By  Remark  2.1,  it  must  therefore  be  a 
viscosity  limit  solution. 

Case  2.  For  Yq  >  consider  the  candidate  limit  solution: 


Uj^(x)  0  <  x  <  1/2 

u_(x)  -  ■ 

u_(x)  1/2  <  X  <  1  . 


Obviously,  we  need  only  check  the  viscosity  inequalities  at  Xq  ■  1/2,  the 
corner  of  U2«  in  this  case  however,  min(u2  ~  can  not  occur  at 
Xq  “  V2  for  any  C^  function  If  max(u2  "  occurs  at  Xq  ■  1/2, 

its  easy  to  check  that  we  must  have  (uj^)jj(1/2)  >  ^^^(1/2)  >  (uj^)jj(  1/2) . 
Computing  these  derivatives,  we  have  that  all  possible  values  of  ^^(1/2)  lie 
in  the  interval  H  ”  Yq.  Yq  ”  1]»  in  which  case 


(t  ‘  »,(4) 


Therefore,  U2  is  a  viscosity  solution. 

Case  3.  For  Yq  <  0,  consider  the  candidate: 

u^Cx)  »  0  . 

Here,  U3  does  not  take  on  its  boundary  condition  at  x  ^  0  or  at  x  1 . 
However,  maxCu^  -  can  not  occur  at  x  =  0  for  any  admissible  test 

,  (-^  (0)  >  “Yg)"  If  on  the  other  hand,  mln(u^  -  is  attained 


function 


at  Xq  ■  0,  we  must  have  that  9,^(0)  lies  in  [YgrO]  and  in  this  case 

(■j  +  UjCO)  >  0  . 

A  similar  argument  shows  that  U3  satisfies  the  viscosity  inequalities  if 
max(u3  -  ^)  is  attained  at  xq  »  1. 

In  these  specific  exaunples,  we  have  demonstrated  that  these  candidate 
limit  solutions  are  viscosity  limit  solutions  of  (3.10)  since  they  satisfy 
Definition  1c.  They  are  furthermore  Lipschitz  continuous  and  so  by 
Proposition  1  and  Remark  3.2  of  Theorem  2,  they  satisfy  ju^  -  u|^  <  const ./e 
where  u^  is  the  exact  solution  of  problem  (3.10).  However,  for  these 
examples,  (as  well  as  other  nonlinear  examples)  there  is  evidence  that 
indicates  a  convergence  rate  faster  than  the  /e,  [3].  We  believe  that  there 

is  a  yet  undiscovered  mechanism  that  links  certain  nonlinearities  in  H  to 
diffusion  vdiich  often  gives  rise  to  a  faster  rate  of  convergence  than 
Theorem  2  predicts. 


54.  NUMERICAL  APPROXIMATIONS. 

In  this  section,  we  Introduce  and  analyze  a  class  of  numerical  schemes 
that  generate  approximations  of  the  viscosity  limit  solution  to  the  one¬ 


dimensional  version  of  (NP^),  which  we  write  here  as: 


du 

-e_+  H(x,u^,—  )  =  0 
dx 


(4.1) 


du. 


dx 


(0) 


'0  ' 


dx 


1  • 


Throughout  this  section,  we  make  the  following  assumptions  concerning 
H(x,u,p),  which  for  ease  of  presentation  only,  is  assumed  smooth. 


Assumption  A*  ;  For  all  xe  [0,1],  |u|  -  R  and  |p|  <  K,  there  exists  a 

>  0  and  an  0  <  <  1,  such  that 

d  ’^1 

H(x,u,p)  >  u,j/(niax(R,  1) ) 

Assumption  B* ;  For  all  x  e  [0,1]  and  |p|  -  K,  there  exists  an 
0  <  t)2  <  1  And  a  constant  C  ( |  u  | )  such  that 

d  ^2 

l-g^  H(x,u,p)|  <  Pj^(max(K,1))  C(|ul)  . 

Assumption  A'  is  merely  a  refined  version  of  Assumption  A  of  Section  2. 

Assumption  B*  guarantees  that  the  viscosity  limit  solution  of  (4.1)  is 

Llpschitz  continuous  and  therefore  supercedes  Assumption  B  of  Section  2. 

The  numerical  approximation  that  are  considered  here  are  built  from  a 

piece-wise  linear  interpolation  of  grid  values  {u.}'^^..  That  is,  we 

J-1  ^  ^ 

partition  the  interval  [0,1]  as  U  [Xj,Xj+^],  where  we  shall  assume  that 

j-0 


2(Xj 


Xj.i)  >  (Xj+,  -  Xj) 


*j-1>  ' 


and  then  define  u^(x)  by 
(4.2) 

where 


u^(x) 


J 

\  UjTj(x)  , 
j-0 


Tj(x)  = 


(x  -  Xj_^ )/(x^  -  X^_^) 

if 

X  e  [X.., 

(Xj+1  -  x)/{x.^^  -  X.) 

if 

X  e  [Xj  ,x 

0 

otherwise 

In  (4.2)  the  superscript  d  is  to  represent  a  measure  of  grid  refinement  and 


we  set  it  equal  to 


max  ( X ^  -  X j ) •  For  each  0  <  j  <  J,  the  grid 
values  Uj  are  required  to  satisfy  the  difference  scheme 

(4.3)  B(Xj,Uj,D'''Uj,D“uj)  »  0 

~D  u  Q  *  Y  Q  *  j  *  Y  ^  t 

where  D'*'uj  =  (Uj  +  i  -  Uj)/(Xj+i  -  Xj),  D“uj  *  (Uj  -  Uj_^)/{Xj  -  Xj_.|)  , 
and  H(x,u,pi,p2)  Is  some  difference  operator  that  does  not  explicitly  depend 
on  any  grid  par^uDeter.  B(x,u,p^,P2)  Is  assumed  to  be  locally  Llpschltz 
continuous  and  It  Is  also  assumed  to  satisfy  three  basic  properties: 

Property  1 .  H(x,u,P),P2)  is  consistent  with  H(x,u,p).  That  Is, 

H(x,u,p,p)  =  H(x,u,p). 

Property  2.  B(x,u,p^,P2)  is  nonincreasing  in  the  p^  argument  and 

nondecreasing  In  the  P2  argument. 

Property  3.  For  all  jp^l  <  K  and  |p2l  <  K,  R(x,u,p^,P2)  satisfies 
Assumption  A'  above* 

IOf  course.  Property  3  simply  says  that  R(x,u,p^,P2)  is  strictly 

increasing  in  u  at  the  rate  prescribed  by  Assumption  A'.  We  now  give: 

t  a 

[  Theorem  3.  With  Assumptions  A'  and  B'  above,  suppose  that  u  comes  frcxn 

[  scheme  (4.3),  where  B(x,u,p‘|,P2)  satisfies  properties  1,  2  and  3.  Then, 

L  L 

(4.3)  generates  a  unique  approximate  solution  u  and  moreover,  u 

converges  to  u  =  lim  u^  at  least  as  fast  as 

c 

I u^  -  u I ^  <  const .  /&  , 

where  above,  u  Is  the  viscosity  limit  solution  of  (4.1)  and 

L  *  max  (’‘j+i  “  X j  ) . 

0^ j<J-1 


Before  proving  Theorem  3,  we  give  two  examples  of  finite  difference 
operators  which  satisfy  properties  1>2  and  3.  Furthermore »  we  show  that  the 
rate  above  is  the  best  possible  under  the  hypotheses  of  Theorem  3. 

Example  1.  The  Lax-Friedrichs  difference  operator  [11],  [20],  is  based  upon 
approximating  h(x,u,^)  by  a  convex  combination  of  H(x,u,d'’’u)  and 
H(x,u,0~u)  along  with  the  introduction  of  an  artificial  numerical  viscosity 
term.  To  be  more  specific,  R  is  given  by 

R(x,u,p^,P2}  -  6h(x,u,p^)  +  (1  -  e)H(x,u,P2)  “  ®^^1  ”  ^2^  ' 

where  6  is  chosen  in  [0,1]  and 

c  >  max(6  sup  IT, (9  ~  1)inf  H  ,0)  . 

P  P 

Clearly,  this  difference  operator  satisfies  properties  1,  2  and  3  ed>ove. 
Moreover,  if  Hp  >  0,  (resp.  Hp  <  0),  we  could  have  chosen  0*0,  (resp 
6  ■  1),  and  c  •  0;  thus  giving  a  scheme  based  on  backward,  (resp. 
forward),  differencing. 


Example  2. 


The  Godunov  difference  operator,  [9] , 

[20], 

is 

given  by 

min 

H(x,u,v) 

if 

P2 

<  Pi 

H(x,u,p^,P2)  ■  - 

ve[p2,p,] 

max 

H(x,u,v) 

if 

Pi 

<  P2  • 

ve[p^,P2J 

This  difference  operator  clearly  satisfies  properties  1  and  2,  and  a 
straightforward  exercise  will  verify  that  it  satisfies  property  3  as  well. 
Again,  when  H(x,u,p)  is  monotone  in  p,  the  scheme  reduces  to  either  a 


backward  or  a  forward  difference  scheme. 


Next,  we  show  that  the  rate  of  convergence  of  Theorem  3  Is  sharp#  We 


again  consider  the  trivial  example  (3.9)  and  we  approximate  its  viscosity 
limit  solution,  (u  «  0),  by  the  Lax-Friedrichs  difference  scheme  -  however, 
we  intentionally  add  too  much  numerical  viscosity,  (we  take  c  =  1  rather 
than  the  allowable  c  *  0).  Setting  -  Xj  “  h,  where  h  =  1/J,  Uj  is 

required  to  satisfy: 

(4.4)  -(D'‘'Uj  -  D^Uj)  +  Uj  =  0 

D“uq  =  0,  D‘*'uj  •  1  . 


One  easily  computes  the  exact  solution  of  (4.4) 

-  ^  f  1  t  ^ 

h  “<k:  “l  ^  -  1  “2  J  ' 


J+1  J+1 

a  -  a 
2  1 


1 


where 


.  h  .  /h  .  (1  .  , 

Oj  .  1  .  I  t  /h  •  (l  , 


and  furthermore,  since  u^  >  0,  we  have  that 


Finally,  calculating  the  right  hand  side  above,  we  arrive  at 


-  “L  >  +  (l  +  ^  =  /h  +  0(h) 


which  is  exactly  the  rate  of  Theorem  3. 


We  shall  prove  Theorem  3  via  three  lemmas. 

Lenma  2.  Assume  that  R(x,u,p^,P2)  satisfies  properties  1,  2  and  3  above. 
Then,  the  difference  scheme  (4.3)  had  a  solution,  say  u^,  u^  is  bounded 
and  has  a  bounded  Lipschitz  constant,  uniformly  in  A  >  0. 


Proof:  We  first  prove  that  must  be  uniformly  bounded.  Suppose  that 

max  u.  >  0  is  attained  for  some  1  <  j.  <  J  -  1.  Since  at  an  interior 
0<j<J  \  ^ 

maximum  u  <  0  <  D  u .  ,  (4.3)  and  Property  2  imply  that 

B(x^  ,Uj  ,0,0)  <  B(x.  ,u^  ,D'*‘u^  ,D“uj  )  ■  0  . 

JQ  Jo  JQ  Jo  Jo  Jo 


Therefore,  we  have  from  Property  3  that 

n.  _ 

U-u .  <  (max(u  ,1))  1h{x  ,0,0,0)]  . 

Jo  ^0  ^0 

Similarly,  if  max  u.  >  0  is  attained  at  j  =  0  or  j  >  J,  we  would  have 
0<j<J  ^ 

that 

- 

yi  [Uq  <  (max(uQ,1))  j  H(  0, 0, 0,-Yjj)  |  , 


or 


yi  lUj  <  (max{Uj,1))  ^|h( 1,0,Y^,0) |  , 


which  proves  that  u  (x)  must  be  bounded  above  independent  of  A  >  0.  An 

identical  argument  would  show  that  min  u.  must  be  bounded  below 

0<j<J  J 

independent  of  A  >  0. 


Next, 

we 

show  that  1 

D'^Ujl 

must  be 

uniformly  ! 

max  0 

f 

u , 

>  niax( 

,0)  is 

attained  at  jn* 

0<j<J-1 

j 

have  that 

0 

a 

'“"“Jo* 

-  H(x^  ,u 

Jo 

> 

'“So- 

D'^u.  )  - 

Jo 

B(x^  ,u^ 
Jo  Jo 

'"jo+1 

'“So' 

-  H(x. 

Jo 

,Ut  ,D'''u^ 

Jo  Jo 

Setting  K  »  D''’u^  >  0,  we  have  from  above  and  Property  3  that 

Jo 


d  ^1 

^  I-5—  H(C,u,  ,K)|(max(|u.  |,1))  , 

^  <»x  Dq  Dp 


and  this  ineqtiality  combined  with  Assumption  B*  implies  that 


^2 

K<  (inax(K,1))  C(  |uj^| )  (inax(  |uj^| ,  1 )) 


Therefore,  D''’uj  is  bounded  above,  again  Independent  of  A  >  0.  A  similar 
argument  would  show  that  D'*’uj  is  bounded  below  independent  of  A  >  0.  This 
proves  the  lemma . 


Lemma  3 .  Assume  that  H(x,u,p^,P2)  satisfies  properties  1,  2  and  3.  Then, 
the  difference  scheme  (4.3)  has  a  unique  solution. 

Proof:  Consider  the  map  defined  by 

(4.5)  (F^(u))^  =  u^  -  vH{Xj,Uj,D‘''uj,D“uj)  , 

for  0  <  j  <  J,  where  “D”uq  =  Yg  and  d'*^Uj  =  We  show  below  that  F^ 

has  a  unique  fixed  point  and  obviously  this  fixed  point  is  the  desired 
solution  of  difference  scheme  (4.3).  We  may  assume  that  I}(x,u,p^,P2)  above 
is  globally  Lipschitz  continuous,  since  R  could  be  modified  in  a  smooth  way 
outside  the  bounded  a  priori  domain  established  by  the  previous  lemma. 

We  now  claim  that  (Fy(u))j  is  a  nondecreasing  function  in  Uj-i»  Uj 
and  provided  that  v  is  chosen  sufficiently  small.  Assume  for 

simplicity  that  H  is  smooth.  We  then  find  upon  differentiating 

(F^(u))  =  VH  /(X  -  X  )  for  1  4  j  <  J  , 

j-1  2 


■5—  (F  (u)).  =  -VH  /(x  -  X.)  for  0  <  j  <  J-1  , 

V  3  P,  3+1  3 

and  Property  2  implies  that  these  quantities  are  nonnegative.  Furthermore, 
(4.6)  I-  (F  (u)) .  =  1  -  v{h  -  /(x.^,  -  x^)  +  H_  /(x^  -  X.,)}  , 


0  for  j  “  J 


for  1  <  j  <  J  -  1,  and  ■  0  for  j  “  0  and 

Therefore,  since  R  Is  assumed  to  be  globally  Lipschitz  continuous,  we  can 

choose  V  small  enough  so  that  these  derivatives  are  nonnegative  as  well. 

Next,  we  show  that  has  the  fixed-point  property  for  v  as  above; 

(V  should  be  thought  of  as  an  artificial  time  parameter  and  the  restriction 
on  V  imposed  in  (4.6)  as  a  CFL  condition).  Let  u  e  and  v  e 

and  define  t  >  v  -  u.  Now  consider 

+  T)  -  F^(u)  . 

Setting  T  »  max[  max  t,o),  we  have  by  the  claim  above,  that 
0<j<J 

(4.8)  (F^(u  +  T)  -  F,(u)),  <  (F,(u  +  T  t)  -  Flu)).  . 

V  V  J  V  W  V  J 

Now,  recalling  the  definition  of  F^  in  (4.5),  we  see  that  the  right  hand 
side  of  (4.8)  is  equal  to 

-  v(H(x.,u,  +  T  ,d\i.,D*u.)  -  H(x  .,u  .,d'^u  .,d”u  . ) )  , 
which  by  Property  3  is  bounded  above  by 

tjjd  '  vu)  , 

A 

where  ^  is  the  appropriate  positive  constant  of  Property  3  governed  by  the  a 

priori  domain  of  Lemma  2.  Setting  t  »  min(  min  T,o}  and  repeating  the 

0<j<J 

argument  above,  we  find  that 

A  A 

(4.9)  T  (1  -  vu)  <  (F^(v)  -  F^(U)).  <  T(1  -  vu)  . 

in  V  V  3  w 

Therefore,  the  Banach  fixed-point  theorem  guarantees  a  unique  fixed  point  of 
F^,  for  V  sufficiently  small,  which  is  the  desired  result. 

Remark  4. 1 .  Inequality  (4.9)  tells  us  that  implementing  an  artificial  time 
method,  (u’^'*’^  =  F  (u*')),  to  obtain  a  solution  of  difference  scheme  (4.3), 


converges  at  an  A  rate  of  e  .  nils  of  course.  Is  computationally  slow 

in  light  of  the  increment  restriction  Imposed  by  (4.6).  We  recommend  a  few 

«• 

iterations  of  artificial  time  to  pull  the  initial  approximation  into  the  I 
domain  of  attraction  for  Newton's  method,  which  with  some  "smoothness", 
converges  at  a  much  faster  quadratic  rate. 

The  next  lemma  is  crucial  to  establish  the  fact  that  u^  satisfies  the 
approximate  viscosity  inequalities. 


Lemma  4.  Suppose  i^x)  =  ic|x  -  yp  +  <|»(x),  where  y  6  [0,1]  is  fixed. 


is  a  constant  and  t|i(x)  is  an  affine  function  with  -i|)'(0)  *  Yq  and 


(1)  =  Y^-  Then: 


(a)  If  K  >  0  cuid  max  (u  (x)  -  ^(x))  is  attained  for  some 

xe[0,1] 

^  0  <  j  <  J,  we  have  D'^’u^  -  D"u^  <  _ (5),  where  x^  is  the 

J  JQ  JQ  XX  JQ 


nearest  grid  point  to  C* 


(b)  If  K  <  0  and  min  (u  (x)  -  ^(x))  is  attained  for  some 

xe[0,1] 

C  /  xj,  0  <  j  <  j,  we  have  D^’uj  -  D“uj  >  A»4>  (C),  where  x^  is  the 

J  Jn  Jn  XX  Jo 


nearest  grid  point  to 


Proof:  We  prove  (a)  only  since  the  proof  of  (b)  is  identical.  Let  x^  be 

Jo 


the  nearest  grid  point  to  C.  We  have  three  basic  cases  to  examine:  They 


are:  x.  =0,  Xj  =1  and  0  <  Xj  <  1. 

30  Jo  ^0 


Case  1:  When  Xj^  =  0,  we  must  have  D'^Uq  =  since  ^  6  (0,x^)  is 


where  the  maximum  of  u  -  4*  occurs.  However,  because  ^  is  quadratic. 


♦^(5)  =  4)^(0)  +  Therefore,  D'*'uo  =  +  (j»^{0)  <  -  Y 


XX 


■^xx 


Case  2:  The  case  when  x.  =1  is  identical  to  case  1  above. 

Jo 


Case  3:  Suppose  now  that  C  6  (x^  i,Xh  )  and  choose  an  arbitrary 

Jo  '  Jo 


T  e  (x^  ,x.  (if  on  the  otherhand,  ^  €  (x^  ,x.  the  argument  below 

Jo  Jq^'  Jo  Jq^^ 


is  essentially  the  same).  Using  the  definition  of  u  and  the  fact  that 


(u  -  ^)(C)  is  maximum,  «re  have  that 


(4.10) 


u.  +  d"u.  (5  -  X.  )  -  4(5) 
Jq  ^0  ^0 


>  u  +  D  u  (t  -  X  )  -  4(t), 


DU.-  4  (5)  . 

JQ  * 

Therefore,  a  simple  calculation  will  show  that  (4.10)  implies 


(4.11) 


4^(5)(5  -  T)  +  4{t)  -  4(6) 

D'*’u  •  ”  D”u4  ^  .  .  ■-  —  — 

30  3o  T  ”  *4 


Taylor's  theorem  allows  us  to  write  the  right  hand  side  of  (4.11)  as 


^0 


Recall  that  we  have  assumed  our  grid  satisfies  the  constraint 

(Xj.|.1  "  Xj)  >  (Xj  -  Xj.^).  Ihis  allows  us  to  minimize  the  bracketed  term 

above  by  choosing  t  -  2Xj  -  5.  Doing  this,  wo  have 


D‘^U4  -  D-u.  <  2(x.  -  5)4_(e)  , 


‘jo  jo  j 


XX 


and  since  x^^  is  the  nearest  grid  point  to  5,  the  proof  is  complete. 


Proof  of  Theorem  3;  The  proof  of  the  theorem  is  conqplete,  (Lemmas  2  and  3), 

A 

except  for  showing  that  u  satisfies  the  approximate  viscosity  inequalities 
of  Definition  2.  With  this  in  mind,  set  4^(x)  -  |x  -  y^l^/i  +  4(x),  where 
yo  6  [0,1],  is  fixed  and  4(x)  is  affine,  with  4'(0)  -  -Yq  and 
♦  '(1)  -  T.|,  (as  in  Definition  2a).  Suppose  now  that  u^  -  4^  is  maximum  at 
5  e  [0,1].  To  show  that  u^  is  an  approximate  viscosity  subsolution,  we 


3 


must  verify  that 


l‘«  It  l‘«  **>  >»»  t*!.  I*,  ISISIVtS 


H{5,u  . 

where  K  is  some  constant.  Independent  of  A,  and  as  always  in  this  section, 
A  =  max  (Xj+1  -  Xj). 

Using  difference  scheme  (4.3)  and  Property  1,  we  have  that  for  every 
0  <  j  <  J 

h(5,u^(5),<fr,„(C))  =  H(C,u^(5),^^j^(C),*^j^{5))  -  B(x^,u^,D■^u^,D■u^)  , 


lj,Uj,U  Uj,u  Uj; 


and  we  rewrite  this  identity  as 


(4.13)  H(C,u“(5),<fr^^(C)) 


=  [H(x^,u^,<|(^^(5)  ,<>^jj(5))  -  H{x^,U^,D  u^,D  u^)] 

+  [H(5,u^(5),t,^(5),^,jj(5))  -  H(x^,u^,*,^(C),*,^(5)))  . 

The  second  term  on  the  right  hand  side  of  (4.13)  is  bounded  above  by 

LxU  -  Xjl  +  '  ’^jl  ' 

where  and  are  the  Lipschitz  constants  of  H  in  the  x  and  u 

arguments  respectively,  and  L  is  the  Lipschitz  constant  of  u^.  The  first 
term  on  the  right  hand  side  of  (4.13)  can  be  written  as 


(4.15) 


where  again,  we  have  assumed  that  H  is  smooth  for  simplicity. 

If  5  =  X.  for  some  0  <  <  J,  we  have  nothing  to  prove  since  it's 

an  easy  exercise  to  determine  that  in  this  case  D  Ua  ^  a  (£)  <  d”u 

-•0  '''lx  ^  jp 

when  X.  ,  (=5),  is  a  maximizer  of  u^  -  (Recall  by  the  definition  of 

Do  ’ 

that  ^  D~Ujj  and  ♦^,^(1)  >  D'''uj  in  the  event  that  5  ”  0  or 

1.)  Therefore,  setting  j  =  jo  (4.15)  and  recalling  Property  2,  (which 


•*  »'■ 


.N 


.*  t''  •** 


Mil 


C'f-Vi 


-2b- 


•ays  that  H  <  0  <  H  )/  varifia*  the  approximate  viacoalty  Inequality  2a 

Pi  P2 

here  in  a  trivial  way. 

If  on  the  other  hand,  ^  ^  for  all  0  <  j  <  J,  take  Xj^  to  be  the 
nearest  grid  point  to  C.  Set  j  ~  jo  in  (4.15)  and  insert  the  identity 
^I^(C)  ■  D~Uj^,  (or  Into  it.  Using  the  result  of  Lenna  4 

allows  us  to  combine  (4.15)  with  (4.13),  to  arrive  at 

H(C,u^(C),«ijj(e))  <  (K*A)*^^j^j^(e)  +  CA  , 

where  K  •  max(-Hp^,Hp^)  and  is  given  by  (Ljj  +  L^L). 

An  identical  argument  will  show  that  u^  is  an  approximate  viscosity 
supersolution,  (see  Definition  2b),  and  so  by  applying  the  abstract  result  of 
Theorem  2,  the  proof  of  Theorem  3  is  conqplete. 
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